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IN FERROMAGNETIC FILMS. 
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Abstract. General phenomenological theory of magnetic spin waves in ferromag- 
netic media is originally reformulated and applied to analysis of magnetostatic waves 
in films and plates with arbitrary anisotropy under arbitrary external field. Exact 
expressions are derived for propagator of linear waves between antennae (inductors) 
and mutual impedances of antennae, and exact unified dispersion equation is obtained 
which describes all types of magnetostatic wave eigen-modes. Characteristic frequen- 
cies (spectra) and some other important properties of main modes are analytically 
considered and graphically illustrated. Besides, several aspects of non-linear excita- 
tions and magnetic-wave chaos are discussed, including two-dimensional "non-linear 
Schrodinger equations" for magnetostatic wave packets and envelope solitons. 



1. Introduction 

This preprint contains some results of research performed in Department of kinetic 
properties of disordered and nonhnear systems of DonPTI NASU, between 2001 and 
2003, and devoted to excitation and synchronization of chaotic magnetostatic spin 
waves (MSW) in ferromagnetic films. In fact, this is a part of report written at that 
time. It reflects our theoretical MSW's considerations which, in our opinion, may be 
useful for concerned readers. Next preprints will reflect our numeric simulations of 
MSW and magnetic-wave chaos in films and corresponding real experiments. 
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It should be noticed that a portion of the following Sections 2 and 3 already was pre- 
sented (in a modified form) by arXiv preprint [Yu. Kuzovlev, N. Mezin, and G. Yarosh, 
cond-mat/0405640 ], but we keep it here for better reader's accommodation. 

2. Basic properties and vocabulary of magnetic waves 
2.1. MAGNETIC TORQUE EQUATION. 

In the classical phenomenology of ferro- and ferri- magnetic solids [1-6] , their magne- 
tization is characterized by 3-dimensional vector field, M(r, t) , representing density 
of elementary magnetic moments. The important property of this vector is that its 
length is constant: \M{r,t)\ = Mg , with = const termed saturation magne- 
tization. Hence, it is natural to measure magnetization and magnetic fields in Mg 
units and introduce the unit-length "spin" vector S{r, t) by means of M = Mg S , 
IS"! = 1 . Then natural time unit will be tq = {2TigMs)^^ , where g is gyromag- 
netic ratio ( g ~ 2.8 MHz/Oe) [1-3]. By this definition, 2'ktq is period of precession 
of a separate spin at field equal to Ms . The rotation of S{r,t) is subject to the 
Landau-Lifshitz torque equation [1-6] which in these units reads 

^ = [F, S] + 7{F - S{S, F)} , F = , (2.1) 

with E being magnetic energy functional and F full effective local magnetic field 
acting on spins (called also local thermodynamic force). Here and below {a,b) , 
[a , b] and a ® 6 denote scalar product, vector product and tensor product of two 
vectors respectively ( c = a ®h is matrix with elements Cq/j = aahp ). After 
transition to dimensionless time, E is factual energy divided by , hence, its di- 
mensionality is spatial volume, while F becomes dimensionless. The term containing 
7 represents "friction", i.e. dissipative interaction between magnetization and other 
microscopic degrees of freedom (thermostat). It can be written in several equivalent 
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forms in view of identities 



F - S{S, F) = [S, [F, S]] ^ {1-S^S)F 



In principle, the Eq.l is the only possible law of magnetic dynamics compatible simul- 
taneously with the constancy l^*! = 1 , energy conservation under spin rotation itself 
and total energy balance with thermostat. In non- principally complicated models 
the friction may depend on spatial gradients of magnetization, VS{r, t) . 

2.2. FERROMAGNETIC ENERGY. 

In any magnet, its energy functional E includes 

(I) anisotropy energy 



which dictates locally preferred magnetization orientation ( dr is volume differential); 

(II) exchange interaction energy, which ensures local spatial uniformity of S{r) 
pattern. 



it is invariant with respect to S rotation; in general, it may be anisotropic with 
respect to gradient; 

(III) dipole interaction energy which causes long-range non-uniformity (demagneti- 
zation) of S{r) pattern and initializes formation of magnetic domains: 




(2.2) 




(2.3) 




(2.4) 



with G being the integral operator. 




(2.5) 
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and G{r) the 3x3 -matrix function 



(IV) interaction with an external magnetic field Hg{r, t) , 




In summary, 

E = E{S, He} = Ea + Ee + Ed + Eh + ... 

The dots replace additional interactions for instance, with distortions and defects of 
atomic lattice. As a rule, their contribution to long magnetic waves related phenomena 
(to be under our interest) can be either neglected or reduced to renormalization of Mg , 
7 or other material parameters [1,2,5]. From experiments it is known that especial 
surface anisotropy may cause pinning of magnetization on the sample surface. This 
effect can be described in terms of boundary conditions for the spin field S{r, t) [8]. 

For example, in good quality yttrium-iron garnet (YIG) ferrites [4,5], Mg ~ 140 Oe, 
correspondingly tq ~ 0.4 ns the exchange interaction radius Tq ~ 5 • 10~^ cm, 
characteristic anisotropy field is less than Mg (that is \dA{S)/dS\ <1 ), and the 
dimensionless friction coefficient 7 is rather small, 7 ~ 10~^ . 

2.3. DIPOLE INTERACTION. 

The dipole contribution (4) is nothing but approximate form of complete magneti- 
zation interaction with self-induced magnetic field, Hs{r,t) . Usually, characteristic 
spatial scales (sample size, wavelength) and group velocity of magnetic waves are much 
less than length and velocity, respectively, of free-space electromagnetic wave with 
the same frequency. Therefore, Hs{r,t) can be obtained as solution to quasi-static 



Maxwell equations 

div (Hs + ^nS) = 0, curl Hs ^ 0, (2.7) 

even in high-frequency (microwave) band. Besides if there are no metal surfaces in the 
vicinity of magnetic sample (waveguide walls or other), one can suppose Hs — > at 
infinity and come to (4). Oppositely, Eqs.7 should be solved under actual boundary con- 
ditions at surroundings and expression (4) be generalized as — — J {Hs S) dr/2 . 

The principal peculiarity of dipole interaction is its specific long-range character: 
with increasing distance |r| it decreases exactly as the volume marked by this distance 
grows G{r) oc . Consequently, the integral dipole force (5) is indifferent to 
absolute spatial scales of magnetization pattern and depends on their ratios only. In 
particular, for uniform magnetization it depends on a shape of magnetic sample but 
not on its size. 

2.4. STATIC MAGNETIZATION. 

Let 5*0 (r) stands for a static (equilibrium, metastable or unstable) magnetization 
configuration in a constant external field, = Hq^t) . According to Eq.l, the essence 
of static state is that everywhere spins are oriented in parallel to thermodynamic force, 
F\\S , that is 

7^ = Wo{r) So , (2.8) 

where scalar function Wo{r) (complete internal magnetic field) is determined by the 
requirement |5'o| = 1 . Generally speaking, this equation possesses many solutions 
(various domain structures) whose physical realization depends on history of the sys- 
tem. But if Hq is sufficiently strong then an unique solution must stay only [1], even 
in spite of the dipole interaction, since the dipole force (5) is bounded, \G S\ < 47r , 
at any nonsingular magnetization distribution. 
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2.5. MAGNETIC EXCITATIONS. 

Consider perturbation of static state, s = S — Sq , for instance, caused by addi- 
tional time- varying external field, h{r,t) , when Hg — Ho+h . Under a perturbation, 
the energy divides into two parts E — E^^ + E± , the ground 

E\i ^ ~J{Wo + {HoSo)}dr - j{h,So)dr, 

and excess energy E± which vanishes at s = and in its turn consists of two parts 

E±{sh} = ^ Jwos'^dr + E{sh} , (2.9) 

E{sh} = J A{s) dr + Eg{s} + Ed{s} - J {h, s) dr 

The first contribution to (9) arises from interaction between s{r,t) and static mag- 
netization (lowering static magnetic order by magnetic excitations when WQ{r) ~ 
(5'oF{5'o}) > , or, in opposite, maintaining it, when Wo(r) < ). The second, 
E{sh} , represents s(r, t) interaction with itself. Its construction is quite similar to 
E{S, H} except that anisotropy term is modified: 

A{s) = A{So + s)- A{So) - {sAf{So)) , Af{S) = 

This self-interaction induces the response operator L , as defined by 

^ SE - 

Ls = — +h = Af{So + s) - Af{So) + {G - rlV^} s (2.10) 

OS 

Evidently, in case of easy axis or easy plane anisotropy, when A{S) and hence A{s) 
are quadratic forms L is purely linear operator. In general, it is reasonable (and 
commonly used) approximation if possible dependence of L on _0 is took into account 
while its nonlinearity (dependence on s ) neglected, thus replacing L by linear 
operator 



L^A- rlV' + G , where (AsU = 441^ " (^.11) 
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2.6. HAMILTONIAN FORMULATION. 

Due to l^o + <s| = 1 , only two of three components of the excitation s{r,t) 
must be considered as independent variables namely, the components perpendicular 
to So{r) , since weak excitation is always perpendicular to So{r) . Thus define the 
two-dimensional vector 

S^^fiS ,U = 1- So^So, 

where matrix H = n(r) performs projection onto the plane, n(r) , perpendicular to 
5*0 (r) . Evidently, 

+ {S\\ -l)So, S\\ = {So, S) , (2.12) 
where scalar represents projection of S onto the static magnetization direction, 
and 

Sl + \S^\^ = 1, |s|2 = 2(1-5||) (2.13) 
In terms of S± the torque equation (1) transforms into 

^ = ^l|[F^5o] + 7(l-^x®^x)F^, (2.14) 
at dD± 

Here Ej_{sh} should be expressed, with the help of (12) and (13), as a functional of 
S± , yielding 

--i(-(^«f))-.flf 

Notice that frictionless version of the Eq.l4 (for 7 = 0) follows from the variational 
principle 



dS± ' 



dv 

+ E^{s, h}] dt (2.16) 



1 + ^1 



It takes canonical Hamiltonian form after the change of variables 
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Si. = QV^-\Q\y^ , S\^ = i- |g|V2 , (2.17) 

where vector Q is also situated in the plane n(r) [6-7]. Simultaneously such a 
choice of variables ensures unambiguous parametrization of S\\ even in case of spin 
flipping when can become negative. 

2.7. DYNAMIC EQUATIONS FOR MAGNETIZATION. 

It is convenient to introduce the rotation (precession) operator, R , and damp 
operator, F , by means of 

RV = [^o V] , fv = f{S^)V = ^n(l- Si_ ® S^)V , (2.18) 
where V is arbitrary vector. Then the Eqs.14-15 takes more pragmatic form: 

^^{R-f{S^)}{{Wo+{Sof))S^-Siif} , f = h-L{S^ + {Sii-l)So}, (2.19) 

to be supplemented by relations (11-13). Note that 11^ = 11 , RU = UR = 
R , R^ — —H , and the latter equality is equivalent to R^ — —1 in planes 11 (r) . 

In generalized version of Eq.l9, direct dipole interaction is replaced, in analogy with 
(7), by self-induced field of the excitation, hs{r,t) , according to the rule hs 
—Gs . In such an approach 

f = h + hs-{A-rlV^)s, (2.20) 
and Eq.l9 should be solved together with magnetostatic equations 



div (/i5 + 47rs ) = , curl hs = , (2.21) 

under classical electromagnetic boundary conditions at the sample surface and sur- 
rounding bodies. 
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2.8. SMALL-AMPLITUDE DYNAMICS. 

If a static state 5*0 is stable, i.e. represents (global or local) energy minimum in 
the space of magnetization configurations then E±{sh} , for small s{r,t) , is posi- 
tive quadratic form. Hence, we have rights to speak about small- amplitude magnetic 
oscillations and waves. Linearization of Eq.l9 results in 

^^{R-^U){WS^-h) ,W = Wo+A-rlV^ + G, (2.22) 
or, in the self-induced field representation. 



^-{R- ^miWo + A- rlV') S^-h- hs} (2.23) 

Let the perturbation and hence linear response vary harmonically: h{r,t) = /i(r) exp(— ia;t) , 
and so on ( i = . Then Eq.23 formally yields 



S± = x{h + hs) X(^) = {i^ + {R- imWo + A- rlV^)}-\R - ^U) (2.24) 

Substituting this expression to Eqs.21, we transform the problem to solving partial 
differential equations 

div /2(/i + /i5 ) = , i2{uj) = 1 + Anxico) , curl hs ^ (2.25) 

Here x A* P^^Y ^ol^ of linear polarizability and susceptibility operators 

[2]. With respect to polarization, x acts as 2 x 2 -matrix performing projection onto 
planes n(r) . It may be thought local spin precession response to given local magnetic 
field, but, from the rigorous mathematical point of view, presence of the Laplasian 
makes it integral (nonlocal) operator. 
2.9. LINEAR FREE WAVES. 
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Omitting in (22) friction and external pump, one can find eigenmodes and eigen- 
frequencies of free magnetic waves (MW): 



Si_ = Ve-^^ , -iuV = RWV (2.26) 

It is sufficient to consider positive frequencies only, a; > , since opposite sign means 

complex conjugating the same mode. Let different modes be enumerated by a set of 
indexes k . Since operator W is positively defined, we can rewrite (26) as 

i^kVk = iW"^RW^'^Vk , Vk = W^'^Vk (2.27) 

Operator in the left equation is Hermitian because iR is Hermitian, hence, the so- 
lutions can be chosen mutually orthogonal and normalized to J {V^Vk)dr — ookSmk 
(with Kroneker symbol on right-hand side). Returning to (27) and (26) shows that 



i J (K; m)dr ^iJ{So [Vk K;])dr = (2.28) 



This gives the orthogonality rule for the eigenwaves. 
The vectors V always can be represented in the form 



V{r) = e^'^(''){a(r) + ib{r)} ,a±b, a±So,b±So, (2.29) 

where both a and b arc real- valued vectors situated in plane 11 (r) and perpen- 
dicular one to another. This means merely that spin vector draws elliptic trajectory 
whose main axes are just a and b . 

2.10. PLANE WAVES AND POLARIZATION DECOMPOSITION. 

Consider free MW in an infinite-size uniformly magnetized sample ( So — const ) . In 
this envisioned situation, eigenmodes (29) are plane waves that is ■j?(r) = {k, r) , with 
k being wave vector, and a and b are constants. One wave only corresponds to 
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every /c , because magnetization precession is definitely clockwise. Putting on Vk{r) — 
Vk exp(i {k, r)) , Vk = + ib^ in the Eq.26, one easy obtains 

- Uk[SQ bk] + iuJk[So at] = Wau + iWbk , W = e-'<'='''> wfte'^'^'^^ (2.30) 

The vectors a and b always can be such ordered that (S'o [a,b]) > , then (30) 
yields 

bk — ■ — ■ l^fcl 
^kPkdk + i^k— = Wak + iWbk , Pk = i — r (2-31) 

Vk \ak\ 

Since operator W , by its definition, is real- valued symmetric operator, this equation 
clearly shows that 

Wak = akttk , Wbk = Pkbk (2.32) 

Hence, the main axes of MW polarization ellipse, and bk , are nothing but eigen- 
vectors of operator (matrix) W . Besides two corresponding eigenvalues determine 
frequency and eccentricity, pk , of MW: 

Pk = V c^k/Pk , ujI = akl3k = det W (2.33) 

(of course, symbol det relates to 2-dimensional projection subspace, while the third 

eigenvalue is zero). 

If the rotation operator R was discarded from Eq.22, instead of —71! , this equa- 
tion would describe monotonic decay instead of oscillations with two relaxation rates 
jak , jPk ■ In this sense, Vk = ak + ibk is decomposition of MW polarization to 
two relaxation modes. 

2.11. SPIN WAVES AND MAGNETOSTATIC WAVES. 

According to (22) and (30), 



W = Ii{Wo + A + + G{k)}li , 
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(2.34) 



G{k) = e-('=''-)Ge^<'=''-) = 47r^ , (2.35) 

where G{k) is Fourier transform of dipole interaction, and A;^ = |A;P . To express 
the frequency more or less evidently, let us introduce the wave vector projection onto 
the plain 11 , k± = Hk , the projected anisotropy matrix A± = UAH and two 
its eigenvalues Ai^2 ■ Then 

ul = {Wo+rie+Ai){Wo+rlk''+A2) + ^ [kl{Wo + r^e + A, + A2) - (^k, UAUk^} 

(2.36) 

ak + ^k^Tr W ,Tr W ^Wo + r^e + + + Airkl/e 

The latter equality as combined with (32) determines the eccentricity. Notice that 
kl^ [Sok]-" . 

Alternatively, the frequency can be found from the Eqs.25, if exclude friction and ex- 
ternal field and represent the magnetization field in potential form, hs (x k exp(i {k, r)) 
This leads to the dispersion equation, 

{k, i2{uj)k) = + 47r (^k, {ioo + R{Wo + A + r^k'^)}-^ Rk'^ = , (2.37) 

whose solution ou — ouk coincides with (36), while wave polarization follows from 
Eqs.24. 

This consideration fully neglects actual size and shape of a ferromagnet sample. 
Dipole interaction in the envisioned boundless MW finds no spatial scales to compare 
with its wavelength. That is why dipole contribution to the dispersion law (36) is 
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insensible to wave length being function of wave direction only. Perhaps the waves 
whose length is ten or more times less than least from the sample dimensions may be 
treated in such a way. This MW are called spin waves (SW). Neither dispersion nor 
polarization of SW depend on sample geometry. 

However, in relatively long wave dipole interaction inevitably compares their length 
with sample dimensions. As the result, dispersion law acquires the ratio of these two 
scales: 

ujk^uj{D\k\,k/\k\,Ty), (2.38) 

where D is some characteristic sample size. Waves which for this correction is essential 
are called magnetostatic waves (MSW). In most of modern apphcations just MSW 
are directly excited and detected while SW generated from them in nonlinear wave 
processes. 
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3. Linear waves in films and plates: propagator, excitation, and 

dispersion equation 

In compactly designed microwave devices current carrying conductors (wires) are 
rather suitable inductors and antennas for magnetic waves (MW) in ferro- and ferri- 
magnetic films. Therefore, it is reasonable to unify consideration of the waves and 
consideration of electromagnetic impedance of conductors interacting with film. We 
will assume characteristic sizes of a circuit be less than 2t: c/ui , with ui being char- 
acteristic operating frequency and c speed of light. This allows to use quasi-static 
Maxwell equations. 

At the same time, real processes in thin films (thin in the relative sense that thickness 
D , is much less than length and width) can be successfully analyzed in the formally 
infinite film approximation [2,3]. In particular, domainless static magnetization of a 
thin film is uniform everywhere except narrow regions adjoining its edges with width 
of order of film thickness. Concretely, evaluation of dipole force (Eq.2.5, that is Eq.5 
in Sec. 2) in finite-size plate geometry shows that demagnetizing field stipulated by film 
edges decreases at least as 2D/d where d is distance from a nearest edge. Thus one 
can surely suppose Sq — const . 

Being relatively thin a film may be thick in the absolute sense that D » tq ( tq 
is exchange radius). We will consider small-amplitude (linear) MW excitation in such 
a film, by wires whose radius much exceeds ro . So thick wires in linear regime can 
directly excite long magnetostatic waves (MSW) only. Therefore, it is reasonable to 
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simplify mathematics by neglecting exchange contribution to the polarizability x (see 
Eq.2.24). This transforms x from integral operator into matrix which locally connects 
spin precession with magnetic field. 

Under these conditions general formula (2.38) for MSW frequency takes the form 



i^k^ i^N{D\k\,k/\k\) , (3.1) 

where k is in-plane wave vector and integer N enumerates wave branches different 
with respect to thickness. Importantly, it shows that (at least at given wave direction) 
the group velocity of MSW, Vg — du^/dk oc , is as small as film thickness. 
3.1. WAVE EXCITATION BY EXTERNAL CURRENTS. 

The external field /i(r, t) (see Sec. 2) produced by currents is solution to the equa- 
tions 

/i(r, i) = ^ hm{r)Im{t) , div /in = , curl hn = —Jnir) , (3.2) 

where , and hnin are total current, its density distribution (normalized to 
unit) and its field for n -th wire, respectively. Clearly, so defined /i„ are nothing but 
Green functions which determine simultaneously electric moving force (EMF) induced 
in conductors by time varying film magnetization. The expression for EMF, > 
follows merely from Eq.2.14 as the consequence of energy balance and reads 

/dS 
{hn -^) dr (3.3) 

This formula is valid for arbitrary strong nonlinear excitation. For weak excitation to 
be considered, EMF is linear function of currents and relation 



(3.4) 
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defines mutual impedances of the wires . 

In linear regime we may treat dipole interaction by means of Eqs.2.23-2.25. Both 
external and magnetization induced field can be represented in potential form, 

/in = — Vt/n , hs = ^ hsmlm hsn = —^Ugn 

excluding (for external field) places occupied by currents. Here Ugn is contribution 
obliged to n -th inductor. The hat means that it is time convolution operator. Suppose 
film situated in region —D/2 < z < D/2 in parallel to xy -plane, and expand 
variables and patterns into Fourier series: 



^n{t) ^ J ^ ^ ' ^n{r) = j e:x.p{ikxX + ikyy)Un{k, z) dk , dk 



(27r)2 ' 

and so on, where k — {kx ky} is in-plane wave vector. The potentials strictly deter- 
mine magnetization: 

S±{uj,k,z) = -x{uj)V^Im{uj)[Usm{i^,k,z) + U m{k,z)] , 

where, after Fourier transform is made, V = {ikx iky V^} . 
3.2. SOURCE FORM-FACTORS. 

In the film interior V^C/n = , therefore external potentials possess simple expo- 
nential behavior: 



Un{Kz)^^n{k)eM\k\{anZ-D/2)}, \k\ = ^kl + kl (3.5) 

where an = ^ (or — 1 ) if n -th inductor is placed above (or under) film, and the 
form-factor, $n(^) , reflects distribution of n -th current. In particular, let n -th 
wire has round cross-section and oriented strictly along y -axis at position x„ , and 
distance between its center line and film top or bottom equals to Pn , then 
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^n{k) = T-j- exp(-|A;|p„ - ikxXn)5{ky) , (3.6) 

with 5{) being Dirac delta- function. 

3.3. IMPEDANCE TO MAGNETIC POTENTIAL RELATION. 

Besides the equality V^C/„ = as combined with usual boundary conditions 
(continuity of normal magnetic inductance and tangential filed) allows to express the 
impedances through the potentials taken at film top or film bottom: 



= e^-*Z„„e-^'^* = — y \k\Usm{co, K a„D/2)t/„(-A;, a„D/2) dk (3.7) 

(we omit details of calculation). Hence, potentials are taken at the surface most close 
to receiving ( n -th) conductor. More exactly, in view of the definitions (3) and 
(4), Eq.8 describes the film contribution to full impedance (which, of course, has also 
contribution from direct magnetic interaction between conductors). 

3.4. WAVE POTENTIAL AND NORMAL WAVE NUMBERS. 

According to Eqs.2.25, we have to solve the equation 

{V,ilV){Us + U) = ,il=l + iTTxioo) , (3.8) 

addressed to either potential of particular current or total one. Of course, the unit here 
(as well as in Eqs.2.21) means nothing but the unit matrix. To carefully consider this 
equation, we need in introducing two 3-dimensional unit vectors 



iy = {k,/\k\,ky/\k\,Q},z = {0, 0,1}, 

and 2x2— matrix 

(z/, /iZ/) (z/, /iz) 
(2;,/iz/) {z'jlz) 



M 



l-^zz 
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The solution to Eq.9 is 

Us + U — u+ exp{X+\k\z) + u- exp{X^\k\z) , (3.9) 
where A± = X±{uj) are two roots of the dispersion equation which follows from Eq.9: 

IJ'zz^'^ + 'i{l^uz + IJ'zu)X- lluiy = , (3.10) 



2 



1/2 



A = A± = Ao±A, Ao = -^^^^^|^^, A^J^-(' ^"^^+^^" ) } (3.11) 

i^zz I i^ZZ \ l^zz 

In general A± can be either real or imaginary or complex. 

Obviously, matrix M and hence A± are sensitive to direction of the in-plane wave 
vector k but indifferent to its absolute value, \k\ . Therefore, as the Eq.9 shows 
potentials and magnetization (5) vary in z -direction as quickly or slowly as in plane. 
This is characteristic property of magnetostatic waves (see Sec. 2). 

To find the coefficients Uj and then boundary values of potentials Usn{i^,k, ±D/2) , 
for substitution into Eq.8, we must once again use standard boundary conditions. The 
final result for the boundary potential, at the side closest to a given inductor, is pre- 
sented by formulas 

Usnioo, k, anD/2) = ^n{k)F{uj, k) , (3.12) 



p. [1- A-z(/i,.-/i.,)]sinh(A|A;|£>) 

^ ' '~ {1 + A) sinh(A|A;|L>) + 2/x^^Acosh(A|A;|L>) ' ^ ' 



A = det M = n^i^Uzz - Huzl^zu (3.14) 

Let us pay attention to the term i{iXzv — IJ-vz) in the nominator of (14). It changes 
sign when k turns in opposite direction, —k . This is example of unreciprocity 
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inherent to phenomena under static magnetization, which leads to difference between 

^nm Qjlld. Zjj-^'fi • 

3.5. POLARIZABILITY MATRIX. 

After excluding exchange from the polarizability x (see Eq.2.24), the inversion of 
denominator operator in Eq.2.24 reduces to algebraic manipulations. The result reads 



{Wo + Ai + A2)U-A^ + iujR . CO 



with Ai^2 being two nonzero eigenvalues of the projected anisotropy matrix A± = 
UAH (see Sec. 2). Additional simphfication is implied by the inequahty 7 << 1 
which is always hoped in practice and, indeed, it is very well satisfied in ferrites ( 
7 = 10~^ -7- 10~^ ). By this reason one can discard all the terms terms containing 7^ 
(all the more higher powers), in particular, make no difference between ou and uj . 
3.6. SUSCEPTIBILITY MATRIX. 

To detail the susceptibility matrix M , a few definitions are necessary. Let 9 be the 
angle between static magnetization vector Sq and z -axis orthogonal to film, and (p 
the angle clockwise counted between y -axis and projection of Sq onto the film plane 
X - y (so that 9 = tt/2 and ip = correspond to in-plane magnetization strictly 
parallel to y -axis). Next, involve to consideration two unit- length eigenvectors of A± , 
Ai^2 (of course, mutually orthogonal), which lie in the plane n_L5'o and correspond 
to the eigenvalues Ai^2 ■ We always can order them so that [Ai A2] = Sq . Then 
consider the plane, z - So , what passes through both z -axis and So , and define one 
more angle, ip , be the angle clockwise counted between this plane and Ai . Thus 
ip = means that vector A2 lies in plane x - y . Also we need in the quantities 



= sin ^ + cos tp = {kx sin tp + ky cos p)/\k\ = cos , (3.16) 
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i/j^ = J^x COS (p — i^y sin (f — {kx cos — kySm. ^)/\k\ = sin , 

which represent cosine and sine, respectively, of the angle, (j) > between the in-plane 
wave vector k and projection of Sq onto plane x - y (i.e. between k and z - Sq 
plane). 

In addition, it is convenient to use the notations 

= (Ai + A2)/2 , A_ = {A2 - Ai)/2 , W^,2 = + ^1,2 , ^1 = W^W^ , (3.17) 



n = Vt{uj) = {Wo + Ai + A2)Ti-A^ + iuR (3.18) 

Recall that at given external field, the vectors Sq , Ai^2 and scalars Wq , Ai^2 are 
completely determined by solution to the static magnetization equation (2.8). In these 
notations Ai — A+ — A_ , A2 — A^ + A_ , and 



where the four matrix elements of Q are defined quite similarly to elements of matrix 
M , 



= , VLyy = {u, Vtvj , Vt^-y = (^z, VLvj , VLy^ = {i', VLzj (3.20) 

After a lot of algebra one can find: 

^zz = {Wo + A+ + A_ cos 2ij) sin^ 9 , (3.21) 
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= [Wq + A+){i/l + i/J cos^ 9) + ^-{(i^ii cos^ 9 - 1^1) cos 2il) - 2iy\\i/± sin 2^' cos 9} , 

(3.22) 

— — iuji'± sin 9 , D,uz — + iuji'± sin 9 , (3.23) 

= {A-iy± sin 21^ - {Wq + A+ + A^ cos 2^') cos 9} sin ^ (3.24) 

Importantly, these complicated expressions always compensate one another so that the 
determinant (15) has the simple pole only, see below. 
3.7. RESPONSE FUNCTION. 

To comfortably express the determinant (15), the roots (12) of dispersion equation 
(11), and the whole function (14) which connects by Eq.l5 the form- factor of current 
distribution, ^n{k) , and magnetic potential of self- induced field, let us introduce the 
characteristic frequencies oUu , oui , UI2 and ui^ , as follows: 

ul = WiW2 + 47r(TFo + A+ + A_ cos 2ip) sin^ 9 , (3.25) 
= W1W2 + 27r(Q,, + Q,,) T 27rV(Q,, + Q,,)^ - {2v^ sin 9YW1W2 , (3.26) 

ujI = W1W2 + 2ti{VL,, + VL^^) + (47ri/xsin 9f /2 = {cul + cuj + (47ri/xsin 9y}/2 (3.27) 
After one more portion of algebra, eventually we obtain: 
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(3.29) 



In this designations the response function (14) takes the form 



F{u;,k) 



(oJq — u;| — 47rc^i/_Lsin 6) sinh [A(a;)|A;|D] 



(3.30) 



{ul - a;2) sinh [A{uj)\k\D] + {ul - uj'^)A{uj) cosh [A{uj)\k\D] 



Notice that all the above defined frequencies ooq , oou , ooi , 002 and 003 , include 
complex factor Wq = Wq — i'juj , hence, they themselves are complex functions of cu 
although with small imaginary parts (and thus weakly depending on ou ). At 7—^0 
all they become real values which characterize spectrum of free MSW (eigenwaves) . In 
particular (see Sec.4), uju is the frequency of spatially uniform spin precession in film. 

3.8. FILM INDUCED IMPEDANCE OF CONDUCTORS. 

As the consequence of Eqs.8 and 13, mutual impedance of two conductors situated 
on one and the same hand from film, is given by 



where, as well as in (8), integration is performed over all the two-dimensional in-plane 
wave vectors while F{u!,k) is determined by Eq. 14 and previous listing of matrix M . 

In all the above formulas the letter uj (as well as all the frequency related designa- 
tions (x!u , i^k , W , W^o , etc.) means dimensionless angular frequency which equals 
to actual dimensional angular frequency 27r/ (or 27r/u , etc.) expressed in units of 
1/to — 27rgMs (see Sec. 2). Hence, correspondence between quantities hke u and 
/ is established by the rule 




(3.31) 



uj = f [GHz]//o[GHz] , /o[GHz] = 5[GHz/kOe]M,[kOe] , 



(3.32) 
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where / is dimensional frequency, g 2.S GHz / kOe is gyromagnetic ratio, Mg is 
saturation magnetization, and square brackets enclose physical unit names. For YIG, 
/o f« 0.39 [ GHz ] . 

3.9. IMPEDANCE OF STRAIGHT WIRES AND LOOPS. 

Substituting (7) in (32), we obtain mutual impedance for a pair of cylindrical wires 
which are directed along y -axis in parallel one to another and to film and lying on 
the same hand from it (again we omit rather tremendous manipulations). To express 
the result in pleasant form, let us measure distances and sizes in centimeters and 
impedances in Ohm. Then 



^nm[Ohm] _ ^ . -q(p„+p^) [(1 - A) cos(gx) + (/i^^ - /i^^) sin(g3;)] sinh(AgL') dq 
w[cm]/[GHz] ^ Jo ^ (1 + A)sinh(AgD) + 2/x^^Acosh(AgD) q 

(3.33) 

Here x — x„ — is distance between the wires and w is their length (in cen- 
timeters). The integral is taken over q = > , ky — , that is in all n -related 
parameters of the integrand it should be put on u = {1,0,0}, z/j_= cos ip , i/y = 
sin (fi , while the frequency ou should be mentioned in accordance with (33). 

For the case when two parallel wires are situated on the opposite parties from the 
film, evaluation of corresponding boundary potentials yields (in the same units): 

^nm ^ 27r^ e-\'i\iPn+Pr.)+i<ix L-MD 2//,,Aexp(Ao|g|L>) 1 dq 

wf 1 (l + A)sinh(A|?|L>) + 2/x,,Acosh(A|g|D)/ |g| 

(3.34) 

(to be accompanied by definitions (12) and (15)). Of course, here u — { sign (g), 0, 0} 
in the integrand. Due to the definition (8), Z = R — 27rifL , where = Re Z 
and L = — Im Z/2Trf play roles of resistance and inductance. 

In fact, what we made is evaluation of impedances per unit length for long line 
inductors with not taking into account disturbance of their parallelity and edge effects. 
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Formulas (34) and (35) can be obviously generalized to strip-shaped wires and to the 
case when transmitting or /and receiving inductor consist of two parallel wires which 
form loop (or many wires which form antennae lattice) and can perform more or less 
wave selection. What is for more complicated configurations to investigate them one 
should return to formulas (13)-(15) and (32). 

By its definition, the distance between any wire and film, p„ , can not be less than 
the wire radius. Closely looking on integrands in (34) and (35) we see, firstly, that only 
such the waves are excited (and contributing to the impedance) whose wavelength 
notably exceeds Pn + Pm ■ Secondly, impedance depends on the ratio S = 2D/ [pn + 
Pjn) but not on D or p„ separately. 

3.10. RESISTANCE OF PARALLEL WIRE. 

Exact analytical integration in (34) is impossible. But in the important special case 
when conductors are lasting along external field, Hq , which lies in the film plane, there 
exists satisfactory analytical approximation. In particular, for magnetic contribution 
to resistance of a single wire (at x — Q ), under neglecting anisotropy, we found the 
estimate 



i?rr[Ohm] 2Ti{2Ti + Ho)fX{f) 



exp {--^ arctanh X(/)} , if /„ < / < MHo+27r) , 



w[cm] [1 - X2(/)] arctanh X{f) 

(3.35) 



X{f) = £—^, U = foVHo {Ho + 47r) , arctanh X = ^ In 

{27^ Jo) 2 1 - A 

In these formulas / is frequency in GHz, Hq expressed in units of Ms , and is 
the uniform spin precession frequency of in-plane magnetized film (with no anisotropy, 
see next Section). Outside of the marked frequency interval, resistance turns into zero. 
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Corresponding wire inductance can be simply estimated at 5 < 0.5 only, and in this 
case 

T IviVG r T T"| 

^^lli^ ^ - (2 + HoM e-^(/) Ei(X(/)) , (3.36) 

w [cmj 

with Ei(^) being exponential integral function [4]. 

Notice that at in-plane magnetization Wq — Hq because in-plane demagnetization 
factors are equal to zero. 

3.11. NUMERICAL EXAMPLES AND DAMON-ESHBACH WAVES. 

In general, the integrals can be obtained by numerically. For this purpose, formula 
(14) is better than (31) because numerical procedure deals directly with matrix (16). 
Some results are illustrated by Figs. la-Id, 2a-2b and 3a-3b. Since we are most inter- 
ested in YIG samples the value Ms ~ 140 Oe was substituted. 

The Fig. la relates to the same case as the Eq.36. To show characteristic influence 
by anisotropy, two pairs of curves are presented, for Ai^2 — and A2 — —Ai — 
Ha — 0.5 , supposing its main axes are x - and z -axes at y -directed field (i.e. 
ip = ). As typically, anisotropy increases frequencies of MSW running across the field. 
Unfortunately, real anisotropy in YIG (function A{S) in Eq.2.2) involves essential 
complications. For the present, in next examples it was neglected at all. Inductors 
lie on one and the same top side of film and last along y -axis while Hq is oriented 
either along y - or x -axis. 

The Fig. lb presents impedance for the loop consisting of two parallel wires which 
continue one another and carry the same current but in opposite directions (see inset in 
Fig. Id). Any inductor parallel to magnetization naturally generates so-called Damon- 
Eshbach magnetostatic waves first discovered in [5]. They run perpendicular to 5*0 . 
Uniquely, their dispersion law can be found (if neglect anisotropy) in simple analytical 
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Fig. 1a. Film obliged impedance (per unit length) of cylindric wire lying on film in parallel to 
in-plane static magnetization, at external field H =3M , anisotropy H =0 or H = 0.5M (left 

US d. d. S 

or right-hand curves), film thickness D=10 )i and wire radius p=20 )i (8= D/p=0.5). 
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Fig. lb. Film obliged impedance of two-element loop with current parallel to static magnetization, 
at the wires separation 1=0.05 cm, film thickness D=10 H =3M , H =0, p=20 
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Fig. 1c. Mutual impedances Z^^ and Z^^ of cylindric wires parallel to in-plane film 
magnetization and distanced by x=0.1 cm from each other, at H|^=3M^, D=10)a., 

p=20)a. . Wire 2 is on the right of wire 1 if look in the field direction. 
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Frequency, MHz 

Fig.ld. Mutual impedance Z^^ of two identical loops parallel to film magnetization 

and separated by x=2l=0.1 cm (with 1=0.05 cm being inter-wire distance in loops), 
at Hp=3M^, D/p=0.5. Thick and thin curves are Re Z^^ and Im Z^^, respectively. 

Loop 2 is onthe right of loop 1 if look in the field direction. 
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form (see Sec.4): 

uJk = ujDE{D\k\) = y^Ho (Ho + 47r) + 47r2{l - exp(-2L>|A;|)} (3.37) 

Maximums of the loop resistance on Fig. lb exactly correspond to Damon- Eshbach (DE) 
waves with lengths 2ml , where / is the width of the loop (i.e. inter-wire distance) 
and m = 1, 3, ... odd integers just as one could expect. Notice that all the spectrum 
of DE waves lies above the uniform precession frequency, Uu — ■\/ Hq (Hq + An) . In 
dimensional form, at Ho/Mg — 3 chosen for these examples /„ ~ 2.68 GHz. 

Figs.lc and Id show mutual impedances of two wires and two loops respectively, 
under the same orientation. Both "from left to right" impedance and "from right to 
left" are presented in Fig.lc. The latter clearly demonstrates violation of the reci- 
procity: every inductor placed above the film top better excites waves going clockwise 
from Ho and 5*0 than inverse waves. To change preferred direction one must remove 
inductors under film. The half of the sum of these two impedances by its magnitude, 
equals approximately to impedance on Fig. la. One may find also that maximums of 
absolute value of impedance in Fig. Id and maximums of resistance in Fig. lb take place 
at the same frequencies and have equal amplitude ratios. In other words loops separa- 
tion strongly influences phase of mutual impedance but slightly its frequency filtering 
characteristics. 

Figs. 2a and 2b relate to wires and loops which are oriented perpendicular to mag- 
netization (see inset in Fig. 2b) and hence excite MSW running along it. The essential 
difference from previous case is that these MSW are irradiated symmetrically (reci- 
procity takes place), and their spectrum lies below a;„ . Since this highest frequency 
responds to least wave number, the group velocity of these waves is negative, i.e. di- 
rected in opposite to phase velocity. 



29 



CD 
O 

CD 
T3 
CD 
Q. 

E- 



400 
300 
200 
100 


-100 
-200 



-300 



1 1 

^— Active, wire 

- Reactive, wire 

Active, loop 
Reactive, loop 


1 1 1 1 1 1 

^ ^ jj^a^^^^^ k V V \ \ 


1 


1 1 1 1 1 1 1 1 


1 



2250 2300 2350 2400 2450 2500 2550 2600 2650 2700 2750 

Frequency, MHz 



Fig. 2a. Impedances of wire and of loop with 1=0.05 cm perpendicular 
to in-plane static magnetization, at HQ=3Mg, D/p=0.5, Mg=140 Oe. 




Frequency, MHz 



Fig. 2b. Mutual impedance of two loops with inter-wire distance 1=0.05 cm separated 
by x=2l=0.1 cm and perpendicular to in-plane magnetization, at H =3M , D/p=0.5. 
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3.12. ROLE OF FILM THICKNESS. 

Let us return to wire parallel to magnetization, and consider Fig. 3a which contains 
the series of resistance via frequency curves for different thickness values. Naively, 
one would predict nearly hnear dependence R{D) if suppose the EMF (voltage) be 
proportional to time-varying magnetization, S± , and to a number of contributing 
spins (that is to thickness), £ oc D\S±\ , while S± proportional to exciting current, 
oc / . But Figs.3a,3b and formula (36) show that, at better excited lower part of 
the spectrum, resistance is almost independent on D (higher frequencies are rejected 
merely by the exponential wire form- factor, see Eq.7). 

What is the matter? The answer comes from Eqs.l and 38: the smaller thickness the 
smaller group velocity of irradiated waves Vg , therefore, the smaller is energy outflow 
from the inductor. Assume this outflow be proportional to product Dvg\S±\'^ , equate 
it to the pumped power RP and combine the resulting relation with £ oc D\S±\ . 
These reasoning yield i? oc D/vg . Taking into account that the ratio Vg/D is 
a function of product the D\k\ or, equivalently, of the frequency only, we get the 
explanation of approximate constancy of resistance. 

From Fig. 3a and Eq.36 it is evident that at D > 1.5p the form-factor becomes 
unimportant and resistance almost independent on D in all the frequency region. 
The strong rise of resistance at frequencies close to fu+ 27r/o reflects fast falling of 
group velocity in this region. Indeed, the Eq.38 implies for DE waves 

Vg = D{{Ho + 27rf -uj^j/co (3.38) 

Hence, the assumption i? oc D/vg implies R oc l/(-f^o + 27r — ///o) what is 
qualitatively confirmed by both Eq.36 and numerical results. This rise becomes better 
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clear if notice that less group velocity means greater density of (excited) states (what 
is highlighted by right-hand peaks condensation in Fig. lb). 

Incidentally, we can conclude that oc I / D , i.e. at given current the more thin 
is film the stronger swing of magnetization and thus the closer nonlinear excitation 
regime. But EMF voltage signal remains approximately the same, even if \S±\ is 
comparable with unit. Therefore, to get greater voltage signal in nonlinear regime, one 
is enforced to make film thicker. 

3.13. ROLE OF FRICTION. 

According to Eqs.l6 the factual dimensionless friction coefficient is not 7 itself but 
product r = 70; (also small quantity). In dimensional form, 

r = ^cu/to = 27r7/ < 27r7/o(f^o + 27r) (3.39) 

The latter estimate relates to in-plane magnetized film where, as we could conclude, 
fo{Ho + 27r) is the upper bound of MSW spectrum (for small anisotropy and long 
MSW as compared with tq ; see Sec.4). 

In infinite-size film under above formal consideration any external source generates 
continuous wave spectrum, regardless of concrete F value, i.e. non-resonant excitation 
takes place. That is why 7 in no way manifests itself in (36). But real finite-area 
film has discrete MSW spectrum. If characteristic frequency separation of excitable 
eigenwaves 5f , essentially exceeds their spectral broadening, 5f » r/27r , then 
it is principally possible to resonantly distinguish them. 

Let us allow that we select only one-dimensional set of DE modes running in x - 
direction and uniform in y -direction ( ky — ). Their separation by wavenumber is 
on order of ~ 27r/d , where d stands for wire length. Hence, frequency separation is 
~ fo{du!k/dk) ■{2n/d) = Vg/d (with Vg being dimensional group velocity). 




Fig. 3a. Film obliged wire resistance, Fie Z^^, via frequency at different 

dimensionless thickness 5=D/p : 0.025, 0.05, 0.1 , 0.25, 0.5, 0.75, 1 , 1 .5, 2, 3.5. 
Wire is parallel to in-plane static magnetization, H =3M (M =140 Oe). 
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Dimensionless thickness, D/p 



Fig. 3b. Dependence of wire resistance on dimensionless film thickness 8= D/p 
at different frequencies (2700:100:3600 MHz) and the same conditions as in Fig. 3a. 
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and the resonance is possible if Fd ^Vg . Expressing Vg from Eq.38, for example, at 
i^o ~ 3 and D ~ 10 [fj] , we obtain Vg ~ Ah'^D/tqUu ~ 10'' [cm/s]. For 7 ~ 
3 • 10~^ and D ~ 10 [/x] , any length d « 10 [cm] occurs sufficiently small! 

In such the case is the theory applicable to real films with d ~ 0.5 [cm]? Yes if 
wave selection, under realistic source form-factor, is not so perfect as was assumed. 
Even if a few modes only with nonzero ky ~ 27r/w are excited in addition to ky 
— , the sufficiently small length easy fails down to 0.1 [cm] or less. Besides in real 
finite-amplitude process directly excited modes transmit their energy to other modes 
by means of non-linear wave interactions. Thus the latters effectively increase friction 
and approach situation to the idealized model. 
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4. Linear waves in films and plates: eigen-modes and dispersion laws 

4.1. IN-PLANE WAVE REPRESENTATION. 

Consider magnetic field created by the magnetization wave 



= Ve-''^' , V = V{z) exp{z{k,x + kyv)} , 



34 



where k — {k^ ky} is in-plane wave vector. For any function G{r) whose 3- 
dimensional Fourier transform is known be G{K) , with K = {k^ ky k^} , and any 
function f{z) , the relations take place as follows 

J G{r-r')f{z')e'^''^^'+''yy'Ur' = e^('=-^+'=-2/) J G{k, z-z')f{z') dz , G{k,z) = J e'^^'G{K)'^ 

If applying this theorem to dipole interaction matrix function defined by (2.5) and 
(2.6), in accordance with (2.35) we have 



G{k,z) = 27r 



{k®k)/\k\ 



kx 


sign(z) 


ky 


25{z] 


- 1^1 



-\k\\z\ 



(4.1) 



where 5{z) is Dirac delta- function. Expression (1) is dipole interaction kernel in 
(A;, z) -representation. In this representation, operator G , the linear dynamical oper- 
ator W defined in Eq.2.22, and the self-induced field take the form 

Gf = / G{k, z - z')f{z')dz' W ^Wo + A + rl{\k\'' -Vl) + G , hs ^ -GV (4.2) 

J -D/2 

4.2. SINGULARITY OF DIPOLE INTERACTION. 

From Eq.l we see that in the long wave limit, when \k\D ^ , all the matrix 
elements of G{k,z) turn into zero except the only singular element, G^z^ 4:7t5{z) , 
which performs purely local connection in z -direction: 



W ^ W^o + ^-roV^ + 47r^O^, ^= {0,0,1} (4.3) 

In this limit, hs{z) =^ —A7rzVz{z) , that is the field vanishes everywhere outside 
film while in its interior in any layer z = const it fully reduces to magnetization of 
that layer. 

This trivial fact of magnetostatics means with respect to spin dynamics that if be- 
ing uniformly magnetized any separate flat layer acts on itself only. Therefore, spin 
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oscillations in different layers can behave independently, resulting in strange property 
of MSW: at \k\D << 1 many MSW branches with different normal wave numbers 
have almost the same frequencies very close to the frequency of uniform spin precession, 

Hence, uju is essentially peculiar degenerated point of MSW spectrum. Although ex- 
change interaction forbids too arbitrary z -distributions and removes exact degeneracy, 
the latter remains important in absolutely thick plates hindering resonant excitation 
of too long MSW. 

4.3. UNIFORM PRECESSION. 

Under uniform precession, all spins in the sample rotate with exactly the same 
phase, that is V = const , and d{r) = const in Eq.2.29. According to (3), in this 
situation, as in case of boundless MW considered in Sec.2.10, operator W transforms 
to algebraic one. As the consequence, clearly, again the polarization decomposition is 
possible. Therefore, we can write 

ujl^detW,W=fi{Wo + A + A7rz®z)U (4.4) 

Since matrix W has exactly the same structure as matrix in Eq.2.34, with 1 in 
place of k/\k\ and zero in place of r^k'^ , in fact the answer is already presented by 
Eq.2.36, namely, 

cul = {Wo + Ai) {Wo + A2) + 4Tr{Wo + Ai sin^ + A2 cos' iP) sin' 9 (4.5) 

The angles 9 and V were introduced in Sec. 3. 6. This expression is equivalent to 
Eq.3.26 obtained when considering MW excitation. 
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The polarization and eccentricity also can be obtained with recipes of Sec. 2. 10 and 
2.11. Of course, in absence of anisotropy the main axes look as a \\z and 6 || [Sqz] . 
In this case, 



p=\b\/\a\ = ^l + —sm'9, {A = Q), (4.6) 

that is polarization ellipse of uniform precession is always in-plane stretched (spins do 
not like piercing the plane). 

For further, let us agree that ^ = and t/j — will mark neglected (or effectively 
negligible) anisotropy or specially oriented anisotropy (see Sec. 3. 6), respectively. 

4.4. DISPERSION EQUATION. 

In fact, the results of Sec. 3 are sufficient to analyze all the variety of free eigenwaves 
in (absolutely thick) film geometry. As usually, their frequencies are nothing but poles 
of the response function, F(a;, k) (see Eqs.3.14, 3.31), at vanishing friction. Hence, to 
get the dispersion equation, we must equate denominator of (3.31) to zero. Elimination 
of friction from formulas of Sec. 3 is achieved merely by returning real quantity Wq 
(static effective field) in place of complex one, Wq introduced in Eq.3.16. Thus we 
come to the dispersion equation as follows 

In view of Eqs.3.22-3.28 and Eq.4, here 0^3 = uj3{k/\k\) is depending on orientation 
of in- plane wave vector while u;^ is constant. 

At given in-plane wave vector, Eq.7 has either one or infinitely many real roots with 
respect to u'^ , determining different MSW branches. For any of roots a;^ = '^Ar(^) : 
we can obtain also corresponding pair of out-plane wave numbers defined in accordance 
with Eqs.3. 10-3.12 and 3.30 : 
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iq±{k) = A±(a;) I A;| = iko{uj) ± A(a;) | A;| , ko{uj) = Xo{uj)\k\/i = ' , (4.8) 

with Qx — ^x{k/\k\) being expressed by Eq.3.15, and ui — ujn (k) . Notice that 
s-nd consequently A;o(^) always are real- valued. 
4.5. TWO TYPES OF WAVES. 

In a bulk wave, by its definition, magnetization harmonically oscillates along z - 
axis i.e. A(a;) takes some imaginary value. Clearly, this is the case if a root of Eqs.7 
belongs to the interval 



< u;^ < ul (4.9) 

Otherwise, A is real and magnetization varies exponentially responding to what is 
usually termed surface wave. However, due to the common scaleless nature of MW 
governed by dipole interaction, characteristic exponents ±A(a;)|A;| , are of order of 
in- plane wave number. At D\k\ >> 1 , such a wave is indeed concentrated in the 
vicinity of film surfaces. But in practically important case, when in-plane wavelength 
is greater than thickness this wave is indistinguishable from bulk one. The Figs. la-Id 
and 3a-3b relate just to such "surface-balk" MSW (see below). 

At arbitrary wave vector k , Eq.7 has infinitely many solutions in the interval (9) 
but no more than one solution outside this interval, i.e. there are infinitely many 
branches of bulk MSW but unique branch of surface MSW. 

4.6. PHASE VELOCITY VECTOR. 

It should be underlined that generally neither bulk nor surface waves are standing 
waves with respect to z -coordinate. The matter is that fix and thus ko{uj) differ 
from zero, and actual wave phase is '&{r) — kxX+ kyy -\-kQZ . In other words 
phase velocity vector in MSW is not in-plane oriented but has also non-zero out-plane 
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component. This fact means that usual attempts to find eigenwaves assuming the 
equaUty g_ = — g+ are wrong. According to Eq.3.25, this equahty takes place at 
special orientation of k or static magnetization only. In particular, for A — , it 
is true if cos ^sin 9 — , i.e. if static magnetization Sq vector is either strictly 
orthogonal or strictly parallel to film plane or if k is strictly parallel to 5*0 projection 
onto this plane. 

4.7. CHARACTERISTIC FREQUENCIES. 

Let us discuss characteristic frequencies defined by Eqs.3.26-28. It can be proved, 
firstly, that the sum flzz+ is always positive and, secondly, the expression under 

square root in Eq.3.27 is always non-negative. Hence, indeed ouf and a;| are positive 
values and both the frequencies Ui and U2 do exist. Thirdly, the inequalities 

uj!<iul<ujl (4.10) 

take place, that is the uniform precession frequency is either immersed into bulk waves 
spectrum or coincides with its edge. This is manifestation of above mentioned fact that 
arbitrary internal layer can undergo autonomous precession. 

Next, notice that uj^ > oui (see Eq.3.28). In view of this inequality as combined 
with (10), the dispersion equation (7) can not be satisfied at cu < cui . At the same 
time, it is easy to see that each frequency from the bulk waves interval (9) can be 
solution to Eq.7 at some appropriate |A;| . Consequently, oui is nothing but lower 
bound of total MSW spectrum. 

It is necessary to remember that all the characteristic frequencies except oUy, only, 
are flowing in the sense of their essential dependence on wave direction, even at A — . 
This is specific anisotropy dictated by flat geometry of dipole interaction. 
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Prom previous reasonings we must conclude that the surface eigenmodes if they 
exist at all, possess frequencies higher than any bulk wave, with lo^ being lower 
bound of surface wave spectrum. But, evidently, the Eq.7 does not have roots at a; > 
max(a;2a;3) . Hence, surface waves do exist at those directions v—'kj\k\ which satisfy 
the condition > uj\ . Then Eq.7 has roots at K{y) < \k\ < oo with lowest 
wave number, k{i') , determined by equality a;| = a;| : 

Dk{u) = {cu',-cuI)/{luI-cuI) (4.11) 

4.8. SURFACE WAVES. 

From previous reasonings we can conclude that the surface eigenmodes if they exist at 
all, must possess eigen-frequencies higher than bulk waves with 002 being lower bound 
of their spectrum. But, evidently, the Eq.7 does not have roots at a; > max(a;2a;3) . 
Consider more carefully the condition for surface modes to exist, cUg > a;| . With 
the help of Eq.3.28, this condition takes the form 

(47rsin 0sin 6)^ > - ^ "^T^Vi^zz + ^ui.)^ - (2 sin 0sin ey{{Wo + A+)^ - 

(4.12) 

( (p was defined in Eqs.3.17). The Eqs.3.22-23 imply 

n,, + n,,^(Wo + A+)(l + sm^(l)sm^e)+ (4.13) 

+A_ [(sin^ (f) sin^ 9 + cos 2(f)) cos 2'0 — cos 9 sin 20 sin 21/;] 

Analyzing these formulas one can see that the requirement (12) can be most easy 
satisfied for nearly in-plane static magnetization Sq and then for waves propagat- 
ing nearly perpendicular to 5*0 . Any surface modes with parallel propagation or in 
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normally magnetized film are clearly forbidden. Prom the other hand, taking exactly 
in-plane 5*0 and = 1 , i.e. strictly perpendicular waves we come from (12) to 
inequality 



27r>|^_|, (4.14) 

which is practically always true. 

Hence, if anisotropy is not extremely strong (in the sense of (14)), then surface waves 
definitely exist in some region surrounding the point 9 — (p — Tr/2 . For given 9 and 
(p in this region, let ouh be their highest (upper) frequency (thus ujh is also upper 
bound of MS W spectrum) . Naturally, Uh is achieved at |A;| ^ oo when cot [A|A;|Z)] 
— > 1 , therefore it follows from Eq.7 that 



= ■ , ■ n\2 = n ■ A ■ n + sm (j)sm9] (4.15) 



(47r sin sin ^)^ \ 2 sin (f) sin 9 
Particularly, for strictly in-plane magnetization and perpendicular propagation, A;_LS'o , 

this expression reduces to 



LJh{9 = 7i/2,(f) = 7i/2) = Wo + A+ + 2n (4.16) 

We can expect this is maximum ouh ■ 
4.9. ISOTROPIC CASE. 

All the algebra becomes much more visual if anisotropy contribution disappears. 
This does not necessarily mean that anisotropy is absent at all. Por example, in the 
easy axis or easy plane case, anisotropy energy (2.2) and matrix A (see Sec. 2, Eq.2.11) 
take the form 



A{S) = Ao{uSf /2, A = Aou(^u, 
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where u is unit vector showing easy (or heavy) axis. Hence, at special static mag- 
netization, when 5*0 II u , the projected anisotropy matrix exactly turns into zero: 
nAU = . 

In absence of anisotropy contribution, the inequality (10) reduces to 



uf = Wl + ^t^Wqv\ sin^ d <ujI^WI\ AtiWq sin^ 9 <ujI = W^ + AttWq , {A = 0) 

(4.17) 

We see that under in-plane magnetization, i.e. at 9 — 7r/2 , (i) 002 — ouu , that is 
all the bulk wave spectrum lies under uniform precession frequency, and (ii) there are 
no bulk waves propagating perpendicular to static magnetization, because the equality 
uji — UJ2 takes place at = sin^ — 1 . 

The condition (12) determining the surface modes region now reads 

sin^ 0sin2 9 > sin^ 9o = , (^ = 0) , (4.18) 

Wq -\- Att 

while their upper frequency is presented by 

9 , V 1 / sin ^0 sin </)sin 6*^^ , , „, 

4 \sm 0sm 9 sm 9o ) 

In accordance with (16), its absolute maximum is max u;^ = Wq -|- 27r . 

The lowest wave number of surface modes is achieved at k strictly perpendicular 

to z - -plane. From Eq.ll we obtain 

mi„C«M = C«(.lS.) = ^-^-^ (4.19) 

Naturally, it turns into infinity at 9 — 9o when surface modes disappear. For in-plane 

magnetization it turns into zero, and then surface waves occupy all the sector 
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It should be added that under in- plane external field ^ = 7r/2 , and the static 
internal field, Wq , trivially reduces to the external one, Wq = Hq , because of zero 
demagnetization. 

4.10. DAMON-ESHBACH WAVES. 

As was mentioned, at ^4 = , in-plane magnetization ( 9 = tt/2 ) and sin^ 
= 1 all three frequencies present in (10) coincide one with another. Hence, in this 
specific case identically A{ui) = 1 , and with accounting for (17) the Eq.7 becomes 
linear equation: 



coth {\k\D) = {ul - - ul) , ul = ul + 8tt^ , col = Ho{Ho + Att) 

Its solution is given by the classical Damon- Eshbach dispersion law (3.38). Otherwise, 
unfortunately, Eq.7 can not be solved in such an evident form. But in wide sense all 
the surface modes can be called Damon-Eshbach waves. 
4.11. MSW CLASSIFICATION. 

In general, dispersion law must be obtained numerically. But also we may treat Eq.7 
be evident expression for \k\ as a function of wave direction and wave frequency: 



sign (a;2 - ujD ^y{uJ^ - ) (a;| - u^) 



\k\D ^Re , ' < ttA^ arctan 

(4.21) 

where N is non-negative integer, the square root should be chosen in upper half- 
plane, while arctan in right-hand half-plane, and the requirement \k\D > serves 
as selection rule for permissible frequencies. Under these conditions N plays no role 
outside of interval (9) while inside it N > enumerates branches of bulk waves with 
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Fig. 4a. Equi-frequency lines of lowest bulk MSW mode (left plot) and surface MSW 
(right-hand plot), for in-plane magnetized film (0=90°), at S.||Y, H =3 and H =A =0 . 

l/z^ a - 

The fat line marks edge of surface wave region, |ky|/|k^|=(47t/HQ) 
The arrows look to higher frequencies. 




Fig. 4b. Equi-frequency lines of lowest bulk MSW mode with frequency less (on the left) and 
greater (dashed line on the right) than co^ (uniform precession frequency) and of surface 

MSW (right-hand plot), for film magnetized at angle 9=45°, with S^ || YZ-plane, 

'Wq=3 and H^=0 . The fat curve marks the edge of surface waves. 
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different normal wavenumbers |/c|£)|A(a;)| . In fact at any A^" two modes can be 
distinguished, witli u) < uj^ and a; > cu^ , to be enumerated as N+ and A^— , 
respectively. 

4.12. MAIN MODES. 

Fig.4a shows equal-frequency lines calculated from Eq.22 for the lowest bulk mode 
0— and surface mode under in-plane magnetization and in absence of anisotropy (for 
details see captures). That are just the two sorts of MSW discussed in Sec. 3 as con- 
tributors to impedances at Fig.2a-2b and la-Id, respectively. 

Wonderfully, their frequency spectra, although lying on opposite hands from uiu , 
at \k\D — > are sewed together in the directions (21), i.e. along the edge of surface 
wave region. Therefore, with respect to sufficiently long waves both the sorts can be 
effectively unified into single main mode. Corresponding expansion of Eq.7 or Eq.22 
gives its dispersion law as follows: 



This equation extends the Damon-Eshbach formula (3.38) to arbitrary propagation 
angles although at small wave numbers. The inequality \k\D <Q.2 is quite sufficient 
to apply Eq.23. 

Fig. 4b shows what does occur if static magnetization is put out from the film plane. 
We see that now main bulk mode 0— becomes strongly separated from surface wave 
but instead the latter well merges with mode 0-|- . 

4.13. EFFECTS OF ANISOTROPY. 

To feel principal influence by anisotropy, let us confine ourselves by special case 
■0 = 0, when one of main axes of 11^411 , namely Ai , lies in z - So -plane while A2 
in X - y -plane. Then 




k\D = (jju + 



nDiAnkl - Hok^) 



(4.22) 



k\y/Ho(Ho + A7r) 
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Fig. 5a. Equi-frequency lines of main bulk MSW mode with frequency less (on the left) and 
greater (dotted line on the right) than co^ and of surface MSW (right-hand plot), for nearly 

in-plane magnetized anisotropic film, at G=85°, H =3 , H =1 .5 and \\i=0 . The arrows 
look to higher frequencies. 




Fig. 5b. Equi-frequency lines of main bulk MSW mode with frequency less (on the left) and greater 
(dotted line on the right) than co^ and of surface MSW (right-hand plot), for nearly in-plane magnetized 

anisotropic film, at 9=88°, H =3 , H =-1 .5 and \|/=0 . The fat curve marks the edge of surface wave region. 
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iul = W2{Wi + 47rsin2 ^) , (V; = 0) , (4.23) 
= WiW2+27r[Wi sin^ 0+1^2(1-008^ ^sin^ 0)]+27r[(iyi sin^ 0-1^2 sin^ 6)'^+ (4.24) 

+1^2^ cos^ e cos^ 0(1 + sin^ 9 - cos^ ^ sin^ 0) + 2W^i W^2 cos^ 9 sin^ cos^ 0] , 
where l^i,2 = W^o + ^1,2 = Wq + ^+ + ^- ■ Instead of (20) one can obtain 

min Dk(u) = Dk(u±So) = — ^ t , ,X = Wi-W2 sin^ ^ , (i/j^O) (4.25) 

47rsin ^ — \X\ 

Clearly, the role of is merely shift of the internal field, and non-trivial effects may 
come from the effective anisotropy field Ha = A_ only. 

For example, consider nearly in-plane magnetized film ( 9 = 85" ) with such an 
anisotropy. Prom Eq.26 it is seen that at > {Wi< W2 ) and 9 close to 
90° the beginning of the surface wave sector remains staying at A; = . That is 
illustrated by Fig. 5a. It shows also that interface between surface mode and 0-|- bulk 
mode is rather sharp, so that the 0— mode (left plot) seems be better continuation of 
the surface sector. Indeed, at ^ — > 90° the picture becomes very similar to Fig.4a. 
Hence, positive anisotropy field. Ha > , (i) produces no qualitative change in relation 
between the two sorts of waves. However, it (ii) allows 0— mode to propagate in x - 
direction, (iii) in accordance with Eq.24, it narrows down the surface wave sector, and 
(iv) rises cuu and thus MSW frequencies (as refiected by Fig. la too). 

In contrast. Fig. 5b demonstrates that negative anisotropy. Ha < {Wi> W2 ), 
(i) causes essential change in relative disposition of bulk and surface modes on A; - 
plane. Namely, the 0-|- mode becomes captured in the bubble immersed into the 
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surface sector (the right-hand edge of this bubble is determined by Eq.26). The matter 
is that at 6' ~ 90" , according to Eqs.24 and 25, 



u1 = ul at sin^ > W2/W1 , cul ^ cul at sin^ < W2/W1 

This results also in (ii) precise angular separation of 0— (left plot) and 0+ modes 
which facilitates to treat, at << 1 , all three modes as single anisotropic mode. 

Besides negative anisotropy (iii) obviously expands the surface sector and (iv) decreases 
uJu and MSW frequencies. 

4.14. LONG- WAVE ASYMPTOTICS AND EXCHANGE CONTRIBUTION. 

In general (at sin^ ^ < 1 ), the dispersion law for long MSW, as compared with 
film thickness can approximated by linear function of \k\D , like in Eq.23. The Eq.7 
yields 



u 



(4.26) 

Hence, MSW modes with greater N are more strongly pressed to the uniform preces- 
sion frequency. This paradoxical fact is the consequence of above discussed singularity 
of dipole interaction. 

Under strictly in-plane magnetization ( sin^ 9 — 1 ) , at N — this expression 
reduces to Eq.23 which unify main modes (may be with the surface one instead of 
0-|- ) . But for N > dispersion becomes quadratic leading to even stronger frequency 
compression: 



uj%^ -col^ ±mDfU - u:l)/{7:Nf , (AT > 0) (4.27) 
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However, at sufficiently large N the exchange interaction enters the game and in- 
creases cjat-i- . At least under in-plane magnetization, the exchange contribution can 
be described by the replacement 

Wo^Wo + rl{\k\'' + 4±) ^Wo + rl{TiN/Df , (4.28) 

with qN± being out-plane wave numbers q± (see Eq.8), for N -th modes in exact 
analogy with how exchange interaction contributes to frequency of boundless wave 
(Sec. 2. 11). In the latter equality, we took into account that for higher order long-wave 
modes ro|A;| is negligibly small while their normal wave numbers q± — ±|A;|A(a;) 
are close to nN/D . 

5. Some aspects of non-linear phenomena and chaos 
5.1. NONLINEAR PROCESSES. 

The unique peculiarity of spin waves (SW) and especially magnetostatic waves 
(MSW) in ferrimagnets (e.g. YIG) is that their relaxation rate, F ~ 5 [ // s ~^ ] 
or even less is very small as compared with other wave excitations in solids in the same 
(microwave) frequency region. So low decay ensures effective generation and nonlinear 
transformations of MSW at small pumping power [1-5] . 

At the same time, usually swing of spin precession remains far from spin flipping, i.e. 
\S±\^ << 1 in Eq.2.19, therefore three-wave and four-wave processes only are of great 
importance. In the firsts either (P) some already excited mode with frequency ouq 
serves as parametric pump for two other modes whose frequencies satisfy the condition 
UJ1 + UJ2 = ooo or, in opposite, (G) two modes mix up one another being the source for 
ujq mode. Here the bracketed letters G and P, abbreviate generation and parametric 
excitation. Among fourth-order processes most important one is the combined G-P- 
process satisfying ujs + 004 = (x!i + (jJ2 ■ 
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If accounting for these processes only, the Eq.2.19 transforms into the approximate 
equation, 

7 Q 

= [So {Wo + L)S^-h\- 711(1^0 + L) S^- (5.1) 
-(^SoLS^-h^ [SoSj] - [SoL{\S^\'' So)]/2- 

-\S^\^[SoLSj_]/2+(^SoL{\S^\^So))[SoSj_]/2, L = l-r^V' + G, 

where three rows contain hnear, quadratic and cubic terms respectively. Higher-order 
terms and all the nonlinear contributions to friction (as well as to anisotropy, see 
Sec. 2. 5) are neglected, and most important entries of h{r,t) are kept only. 

Clearly, the external field, /i(r, t) , also can act as either additive source (G-process) 
or parametric pump (P-process). The first variant is more effective if realizes by way of 
ferromagnetic resonance (FMR). In best real YIG samples the power consumption of 
order of tens microwatt may be sufficient to initiate nonlinear processes [2] . At greater 
pump, one can observe rich variety of nonlinear phenomena including formation of 
envelope solitons [2,3-9], parametric amplification [10-12], magnetization reversal [13], 
self-focusing of MW beams[14], generation of harmonics subharmonics and ultra-short 
pulses [15], non-hnear short electromagnetic waves [16] (the alternate to MSW high- 
frequency branch of mutual magnetization and EM-field hybridization). 

But most interesting phenomenon is magnetic chaos (chaotic oscillations of magne- 
tization pattern) produced if external pump exceeds certain critical level [2,17-19]. 

5.2. NONLINEAR WAVES. 

In special class of nonlinear phenomena qualified as weakly nonlinear magnetic waves 
a narrow region of total MSW frequency band (all the more of whole MW spectrum) 
is involved only and, hence, third-order processes (quadratic terms in second row of 
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Eq.l) are not at business. In sufficiently long waves their dispersion (spatial derivatives) 
also is weak and therefore naturally separates from nonlinearity, so that approximate 
wave equation turns into sum of spatially non-local (differential) linear terms and local 
nonlinear (cubic) terms [3]. 

If speak about films the long magnetostatic nonlinear waves are of special interest 
composed by the main branch of linear MSW (most homogeneous with respect to 
normal 2; -coordinate) . In this case the exchange part of operator L can be neglected, 
and dispersion is completely determined by dipole interaction (see Sec. 4). But since 
the singular part of dipole interaction (Sec.4.1) is factually local, its product with cubic 
nohnearity should be kept. Then the Eq.l (as combined with Eq.4.1) reduces to 

^ = [^o (1^0 + ^ + G) S^-h]- 7n(W^o + A + A7rz®z) (5.2) 
at 

+^\S^\^[So(A\\+A7rSl-A-A7rz<»z)S,], A\\ = (^SoASo) , 

where characteristic frequency of h is supposed the same as carrying frequency of 
S± . Of course, still this is formal storage only for more correct equation which must 
be free of third-order harmonics and concern S± 's envelope. Such the equation can 
be deduced, as usually [20], from variational formulation of Eq.l, or by means of time 
averaging over the carrier period. 

5.3. NONLINEAR WAVE EQUATION. 

In accordance with Sec.2.9 and Sec.4.3, the main (as well as any other) branch of 
eigenwave modes looks as 



S± cx Vfc(r) exp{-iLJo{kD)t} , Vfe(r) = {ak{z) + ibk{z)} exp{i {k, p)} 
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Here its dispersion law is written in the form uj — uJo{kD) , k — {k^ ky} is 
in-plane wave vector, p = {x,y} , and ak{z) and bk{z) are mutually orthogonal 
real-valued vectors. Arbitrary non- autonomous (externally influenced) wave composed 
by these modes can be expanded into Fourier integral 

S±^Re J {ak{z) + ibkiz)}e'^'''P^C{k)^{k, t) dk , (5.3) 

where function \l/(/c,t) contains one-signed (e.g. positive) frequencies only, that is 
represents an analytical signal. The C{k) in (3) being real positive factor serves for 
suitable normalization of the eigenmodes. If it is fixed then, instead of (3), one can 
equivalently consider the "wave function" 

* = ^(x, y, t) = e'^'^'f'^^ik, t)dk (5.4) 

Correspondingly to (3), it useful to introduce analytical signal, h , for the external 
pump too: 

/i(r, t) = Re h{r, t) , j e^*h{r, t)dt = at u < 

For the wave function, the Eq.2 implies the equation as follows (we omit its deriva- 
tion) : 

— + iujQ(-iDVW = -i>i\'^\'^'^ -T"^ + r) (5.5) 
at 

Here V = {d/dx,d/dy} and operator a;o(— i-DV) (formally differential) is deter- 
mined by the dispersion law. Let a, b and a, /3 be the pair of eigenvectors and 
related eigenvalues of the uniform precession operator, W= Il{Wo+A+4:7rz®z)U. , 
considered in Sec.4.3, and p eccentricity of uniform precession. Besides for any two 

vectors u and v , let Uy means u 's projection onto v , i.e. Uy = {v,u) /\v\ . 
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Then the parameters of Eq.5, friction coefficient F , nonhnearity scale x and pump, 
1] , read 



r=(a + ^)7/2, (5.6) 



f [VP^a-'^^b/Vp} , P^^, (5.7) 



2{a + P)\, (5.8) 



while approximate connection between the wave function and magnetization is estab- 
lished by 

'^ = VpS±a-^S±b/VP, M'-\\S±a\\\\S±b\\ (5.9) 
Here ||...|| denotes envelope (amplitude) of an oscillating variable. 

In particular case of tangential magnetization and not strong anisotropy, vector a is 
nearly parallel to normal z -axis vector b lies in the film plane, and formulas (6-8) 
are simplified to 

r^{Ho + 27r)7 ,p^y/l + iir/Ho , x ^ - ""[f ° f ^ (5-10) 

Ho + 47r 

5.4. NONLINEAR SHRODINGER EQUATION. 

For waves and wave packets formed by a narrow set of in-plane wavevectors concen- 
trated about some ko , the Eq.5 reduces to the nonhnear Shrodinger equation (NLS), 



+ (vgV) ijj = i (V, DV^ ijj - iH\'il]\^'il^ - + ^, (5.11) 



tp = exp{iuJo{Dko)t - i {ko p)}^' , rj = exp{iuJo{Dko)t - i {ko p)}r] , 
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^ dko ' 2 dkomdkon 

Here Vg and D are group velocity vector and diffusivity tensor, respectively, and ip 

plays the role of envelope of wave function ^ . 

Evident analytical expressions for the latter quantities can be obtained in a few spe- 
cial cases only, particularly, for Damon- Eshbach waves (see Sec. 3 and Sec. 4) in exactly 
in-plane magnetized film with zero (or weak) anisotropy. In this case, if magnetizing 
field Ho is oriented along y -axis then for surface waves nearly parallel to x -axis 
{ ky <^ \k\'^ ) the Eq.4.7 yields (in the dimensionless time units): 

oj^iDk) ^ujI + 7r{l - exp(-2L>|A;|)}(47rA;^ - Hokl)/\kf (5.13) 

Hence, wave packet running along x -axis with koy = , has group velocity and 
diffusivity as follow: 

Vgy^O, Vgx = 47r^ D exp ( " 2i:> | /cq | ) slgn ( /cq) /^Q , <^o = ojDE{D\ko\) , (5.14) 
with function u!DE{D\k\) given by Eq.3.38, and 

3^^ = -4:TT'^D^{{Ho + 2Tif -2Ti\xp{-2D\ko\)}exp{-2D\ko\)/uj^, (5.15) 

D,y^O, Dyy = -7r(//o + ) { 1 - eM-'^D\ko\)}/2\ko\^ (5.16) 

In contrary to this specific case, generally propagation direction of envelope of the 
wave packet differs from its carrier wave direction, ko . Clearly, the group velocity 
is perpendicular to equi- frequency curves shown at Fig.4a-b and Fig.5a-b (see Sec.4). 
These figures (as well as formulas of Sec.4) show that wave packets which are formed 



54 

by surface MSW and have non-zero koy comparable with ko^ must prefer directions 

characterized by 



Vgy p:i ±Vg^^/H^J4^^ (5.17) 

If the carrier wave is not long, that is the value D\ko\ is comparable with unit, then 
the main axes of polarization ellipse, a, b , its eccentricity, p , and the eigenvalues 
a, (3 in Eqs.6-10 should be calculated just for the /cq mode (instead of uniform one), 
i.e. mentioned as Ofeg , hk^ , and so on, in the sense of Sec. 2. 9-10. 

5.5. NON-ISOCHRONITY AND INSTABILITY OF MAGNETIC WAVES. 

Consider autonomous waves i.e. in absence of pump and dissipation. Cubic nonlinear 
terms in Eq.5 and Eq.ll involve fundamental non-isochronity property of nonlinear 
MW: their frequencies depend on their amplitudes. Indeed, for a plane autonomous 
wave with amplitude A the Eq.5 gives 

* = Ae-x.Y>{-i[oJo{Dk) + xA^]t + i {k, p)} (5.18) 

According to Eq.lO, in tangentially magnetized film intensification of wave leads to 
lowering its frequency. 

What does occur if the amplitude is not uniform but slightly spatially modulated? 
As in general [20], result depends on concurrence between nonlinearity and dispersion 
which in oure case is described by diffusional term in Eq.ll. To see the result, let us 
search for evolution of the wave envelope in the form 

i'ip: t) = [^ + X{p - Vgt, t)] exp(-i>f A^t) , x = Xi + 'i'X2: A = const , (5. 19) 

with X being (infinitely) small non-uniform perturbation. It is easy to derive from 
Eq.ll the linearized equations for Xi(p, t) and X2{p,t) as follows 
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instance, Xi(p, 0) oc cos (A;,p) , X2(p, 0) = . Then solution to (20) consists of two 
definitely weighted exponents exp(±Ai) , where ±A are eigenvalues of the right-hand 
matrix operator, 



Here (^q, DqJ represents deviation of wave frequency coming from the dispersion. 
If it is of the same sign as the deviation due to non-isochronity then A has imaginary 
value. Hence, in this case initial non-uniformity results in small amplitude and phase 
oscillations which in reality decay due to dissipation. But if 



then the sufficiently intensive wave occurs unstable with respect to small amplitude 
disturbance. The latter grows and the wave inevitably breaks into a chain of energy 
slots (solitons). 

According to the instability conditions (22), if friction was absent then sufficiently 
smooth spatial modulation always would unstable. Due to friction, however, both long 
and short modulations always are stable (taking into account that in reality A <1 , 
because of relations (9)). The instability starts from moderate modulation scales qc , 
and after exceeding at least minimum threshold amplitude value, ^min , as determined 




(5.21) 





(5.22) 



by (22), 




(5.23) 
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5.6. SOLITONS. 

Magnetic envelope soliton is a single spatially local wave packet stabilized (protected 
from diffusional bleed) by non-linearity and described by Eq.5 or Eq.ll. Its existence 
is implied by the same instability (they must satisfy the first of the inequalities (22)), 
but to created it one should use suitably localized external pump, instead of a spreaded 
wave. These solitons are called also "bright sohtons" . 

The envelope of one-dimensional (flat) autonomous (at no pump at no friction) 
bright soliton, moving in direction of some unit-length vector n , is determined by the 
equations 

^{p, t) = e-^^*r(0 ,^ = {np- Vgt) , 5^ -HF^ + nr^O, S=(nDn) , (5.24) 
which directly follow from Eq.ll. The solution to (24) is 

r{C) = A/ cosh(A{|x/25|^/2) , Q = KbA^jl , (5.25) 

with the magnitude A being free parameter. 

Alternatively, so-called "black" (dark) envelope sohtons can exist representing "holes" 
(dips) in amplitude of spreaded (plane) wave. From the point of view of above con- 
sideration, these objects formally correspond to imaginary modulation wave vector, 
g — )■ ig . Thus for them the first of the instability conditions (22) turns into oppo- 
site, but the second remains valid. Hence, sufficiently intensive MW inevitably loses 
stability and produces some soliton structure by either one or another way. 

The envelope of black soliton satisfies the same equations (24), but with non-zero 
boundary values at infinity, and has the form 



r(0 = Atanh(A{|x/25|i/2) , ^ = kJ^ 



(5.26) 
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Clearly, both the types of sohtons are as much narrow (wide) as strong (weak). The 
peculiarity of the black soliton is phase slip in its center by tt . For a given direction 
n , either bright (if Sx < ) or black (if Sx > ) solitons exist only. In reality, 
two-dimensional solitons are under use [3-9] , but their analytical investigation is much 
more hard task. 

5.7. MAGNETIC CHAOS. 

The comprehend reviews of experimental data on magnetic chaos its theoretical 
interpretation and numerical reproduction are presented in [2,17,18,19]. The important 
conclusion from both theory and numerical simulations is that even two MW modes 
(i.e. four variables: two amplitudes and two phases) are sufficient to realize chaotic 
behavior. 

In standard scenario, the mechanism of chaos is dependence of frequencies of the 
modes on amplitudes because of their nonlinear self-interaction and mutual parametric 
interaction. At small amplitudes they are coherently (resonantly or parametrically) 
excited by external field and one by another. At large amplitudes the coherence de- 
stroys and dissipative damping prevails which restores coherent interaction and returns 
to beginning of the cycle. Under sufficiently strong pump, this cycle becomes unstable 
with respect to infinitely small perturbation and thus chaotic. 

Hence, the same property (non-isochronity) of magnetization oscillations at the same 
degree of nonlinearity is responsible for both chaos in only two-mode model and for 
regular soliton structures consisting of very many MW modes (let us recollect that 
autonomous NSE dynamics is integrable and thus can not produce chaos [20]). This 
fact demonstrates that principal origin of magnetic chaos is nothing but energy transfer 
through magnetic system (from external source to thermostat). In other words this 
is dissipative chaos characterized by phase volume contraction and dissipative strange 
(zero Lebesgue measure) attractors [21], although (due to small friction) possessing 
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many properties of Hamiltonian chaos [21]. Then, it is not surprising that just the 
power absorption (energy consumption by a ferromagnet sample per unit time) mostly 
highlights magnetic chaos [2,17,18] and can be used as control variable for its identifi- 
cation and synchronization [19]. 

In principle, in presence of periodic perturbation accompanied by dissipation even an 
individual spin (magnetic moment with fixed length and thus two independent variables 
only) can undergo chaotic behavior [22-24]. Nevertheless the two MW modes are too 
few to adequately imitate real magnetic chaos since they are forced to incur roles of 
other modes. Therefore a variety of many-mode models was suggested for numerical 
investigation [2] which are able to reproduce (i) typically observed spikes in the power 
time series (ii) their intermittency, (iii) their fractal properties and, moreover, (iv) 
characteristic frequencies of chaotic power oscillations usually in the interval from 0.1 
MHz to 10 MHz. 

5.8. FRACTAL DIMENSION AND CONTROL OF CHAOS. 

Physically, the peculiarity of chaos (in contrary to noise) is that very many degrees 
of freedom are governed by a few independent variables only. What anybody needs 
in when describing chaos is adequate choice of such the relevant variables (which may 
differ from some particular modes). At least the number of relevant variables drei , 
can be determined if estimate so-called fractal dimension, dfrac , of time series under 
observation. 

The quantity dfrac characterizes dimension of a manifold (attractor) filled by tra- 
jectories of the relevant variables. It is obvious that dfrac < drei ■ At the same 
time, dfrac > drei — 1 5 since the opposite case would mean that one of variables 
is somehow dependent on others. For example, if the attractor was periodic (limit) 
cycle whose dimension dfrac = 1 it would be described by single variable (its phase). 
For chaotic (strange) attractor, its dimensionality dfrac is inevitably non- integer. 
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This means that its intersection with (almost any) one-dimensional line (in drei - 
dimensional embedding space) represents so-called Cantor set. The latter is infinitely 
rarefied (nowhere dense) set of uncountably many points [21]. Roughly speaking, if 
line contains points then a Cantor set on it contains 'i^^ points with S < 1 . Then 

dfrac = drel — 1 + S . 

If dfrac is known then the number, drei , of variables which are governing chaotic 
dynamics can be found as the integer number exceeding dfrac but most close to it. 
The cases when dfrac > 3 (and thus drei > 4 ) are called hypcrchaos. For the same 
purpose of drei determination, the correlation dimension can be used. 

The correlation dimension [25], dcor , characterizes statistics of distances between 
points of the attractor taken at discrete time moments tn = + nr , with some 
reasonable time interval r and n = l..A^,at N ^ oo . Let Xd{t) , d — l-.drei , 
be relevant variables under consideration. Then the set of distancies 

d 

is investigated as characterized by the so-called correlation sum. 



Due to fractal (scale-invariant) structure of Cantor sets it can be expected that at 
small distances and large number of points the power law takes place: 



(J^ \ dcor ^ 
— j at N ^ oo , — > , /2max = max Rij , (5.27) 

where Q is some constant, and the right-hand limit of the correlation sum presents 
definition of dear - Naturally, under rather general assumptions dcor = dfrac [25-28]. 



60 

There are two important statements. First, dcor is insensible to smooth transfor- 
mations of attractor variables including (not too long) time delays. Therefore, equiv- 
alently one may analyze discrete sequence, x{to + nr) , of any available variable, 
x{t) (of course, well connected to attractor), considering the subsequences {x{tQ-\-nT), 
x{to + fiT + t), ...,x{to + nr + dr)} , with d > drei , quite like the attractor points 
above. Second, in principle, estimate of dcor is insensible to d (called embedding 
dimension) if only d > d^ei ■ But it is sensible to noise, either external parasitic one 
(errors of measurements etc.) or ham noise produced by a real system itself. Hence, a 
factual dear 's dependence on d can inform about quality of data under analysis. 

Prom the other hand, in general structure of chaos (strange attractor) may be better 
characterized by a spectrum of fractal dimensions instead of a single one [27,28]. Then 
different variables may give more or less different correlation dimensions. 

The essence of chaotic motion is its exponential instability, that is exponential growth 
of response to arbitrarily small disturbance. Nevertheless this motion obey determin- 
istic law. Hence, if its current state is controlled with accuracy up to n binary digits 
then its future can be somehow predicted for a time, n ■ t-ud , proportional to n . 
Then hch = ln2/tinf is called entropy of chaotic attractor, while the sense of tjnf 
is lifetime of information bits. The latter approximately coincides with characteristic 
correlation time of chaotic variables [21,33]. 

Let there arc two identical chaotic generators initially delivered in the same state 
to some extent of precision. To keep the same equality of states in future and thus 
synchronize one generator by another, we should send from one to another at least one 
bit of information per time tjnf . A representative chaotic variable carry just such the 
amount of information and hence can be used for the synchronization [29,30]. The real 
example of magnetic chaos synchronization was reported in [19]. 
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However, to make this minimum necessary information to be also factually sufficient, 
it should be chosen and applied in adequate way. Concretely, one must take into 
account the topology of attractors i.e. graph of transitions between its Cantor subsets. 
For example [34] , the discrete-time chaotic evolution described by the tent map, x{t + 
1) = 1 — \2x(t) — 1| , with < a; < 1 , has hch = ln2 , i.e. one bit of information 
per time step is sufficient for synchronization. But this principal possibility turns into 
reality if only the bit is chosen be at x < Xq and 1 at x > Xq , with certainly 
Xq = 1/2 . Any other rule (or other xq ) cither leads to errors or requires additional 
information. Generally, determination of adequate rule (termed generating partition 
of phase space) and corresponding most meaningful information sequences (so-called 
symbolic dynamics) is very non-trivial task [34,35], even if dynamic law of chaos is 
known, all the more if it is under question. From this point of view, the results of [19] 
seem extremely interesting. 

In principle, an adequate rule allows to synchronize non-identical attractors too 
if they have similar topologies and equal entropies. At more simplified approaches to 
synchronization (but instead practically applied ones see [31,32] and references therein), 
rather small non-identity of "master" and "slave'' chaotic systems can forbid it, even 
in spite of quantitative excess of information. 

Let X(t) and Y{t) are vector ( drei -dimensional) variables of two chaotic systems 
which obey the same dynamic equations but the second is influenced by the flrst as 
follows: 

dX/dt = F(X) , dY/dt = F{Y) -g-(Y-X) (5.28) 

Here g is positive matrix, hence, it introduces additional damping. Let the latter is 
so strong that suppresses exponential instability of the slave system. Then it easy to 
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see that after some time the only possible solution for Y(t) will exactly reproduce 
X{t) , and thus one can say that the slave system is ideally synchronized by the master 
system. 

Perhaps however, it would be more correct to name this copying chaos. Indeed, the 
slave factually loses its autonomy (since at X = it would produce neither chaos nor 
any other motion instead tending to a stable state), and behaves as passive repeater 
of external signal. 

In more general and fine variant of such kind of synchronization, 

dX/dt = F{Xi X) , dY/dt = F{Xi Y) , (5.29) 

where Xi is some (say first) of drei attractor variables and again function F{Xi Y) 
of two arguments is arranged so that first equation produces chaos while solution of 
the second falls into stable point solution as Xi = . 

The copying of chaos is rather sensible to non-identity of the slave and master systems 
to adding external noise or any distortion of master signal in transmission channel (see, 
for instance, [36] and references therein). Relative error of the reproduction occurs be 
at least the same as relative difference of master and slave parameters plus noise to 
signal ratio and plus relative distortions. 

In practical applications more reasonable approach may be to surely recognize and 
reproduce some particular characteristics of chaotic signal only, instead of its literal but 
erroneous copying. The example is mutual phase synchronization of chaotic oscillators 
(for instance, famous Rosslcr systems) which does not need in simultaneous amplitude 
synchronization and therefore is possible for non-identical oscillators in presence of 
noise. More general possibility is so-called event synchronization where events mean 
definite well characterizable fragments of chaotic trajectory. 
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Then the natural step is an artificial creation of events in master system which 
can serve for encoding information and then its decoding in similar slave system. In 
particular, this may be switching between different trajectories on the same attractor. 

When artificially manipulating trajectory of a chaotic system, one needs in a set 
(alphabet) of easy creatable and identifiable "events" . Such the possibility is ensured 
by unstable periodic orbits (UPO), i.e. periodic trajectories which always take place 
on chaotic attractors. Moreover, from practical point of view one may treat a strange 
attractor merely as a collections of periodic orbits with different length, from some 
minimum period up to infinity. The instability of finite-length orbits means that their 
measure (relative number of attractor points belonging them) is zero, therefore, almost 
even insignificant deviation from short periodic orbit for certain injects to very long 
one (chaotic). 

But, remarkably, UPO's can be stabihzed and thus practically installed into master's 
chaotic trajectory by means of specially programmed feedback (see, for example, [37]). 
Then similar feedback in slave system helps to unambiguously recognize an UPO's 
installation although it looks quite as typical fragment of transmitted signal. Such the 
discrete chaotic encryption of information can well protect it from noise and signal 
distortions. 

During recent decade many ideas of chaos application to secure communication were 
suggested. One of schemes successfully realized in [31] is based on introducing com- 
munication signal, s{t) , into Eqs.29 : 

dX/dt = F(Xi + s{t),X) , dY/dt = F(Xi + s(i) , Y) (5.30) 

Thus the master (in [31] it is chaotic Chua's generator) produces chaos influenced by the 
signal. The factually transmitted information is the sum Stransif) = ^i{i) + s{t) . If 
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the communication s{t) was absent Y{t) would be precise copy of X{t) . Therefore 
the communication can be restored as s(t) = Stransif) — Yi{t) . 

Analogously, discrete chaos (chaotic maps) can be used. For simplest example, let 
x{t-\-l) — F{x{t)) be some one-dimensional chaotic map, and we introduce discrete- 
time information s{t) by means of x{t + 1) = F{x(t) + s(t)) . If Stransit) = 
x{t) + s{t) is sent, then in identical receiving system the information can be recovered 
merely as s{t) — Strans{t)— F'{strans{t — 1)) ■ This possibility was suggested in [38] 
for information encoding in chaotic impulse communication. 

Principally similar ideas were experimentally realized for communication with optical 
chaos [39] . The peculiarity of the latter is an essential time delay in a feedback part of 
optical (laser) chaotic generators. Correspondingly, their dynamics undergo difference- 
differential nonlinear equations which can produce chaos whose fractal dimensionality, 
dfrac , exceeds formal number of variables (number of equations). 

Such the schemes (in which an information is either masked by chaos or modulates it) 
possess all the potential defects of chaos copying. Besides they do not allow for multiple- 
access chaotic communication (many senders) in the same time-frequency domain. The 
matter is that chaotic system can not recognize even its own signal if is mixed with a 
signal from other system. 

Still chaotic extension of modern digital code-division multiple access (CDMA) com- 
munication is under discussion. Although CDMA also uses chaotic signals (pseudo- 
random coding sequences) but these are discrete exactly predictable (periodic) signals 
only whose entropy is zero. The interesting scheme of multiplexed chaotic commu- 
nication based on analog chaotic signals was suggested in [40,41]. It shows how the 
unrecognizability of mixed chaotic messages can be overcame. The idea is that all 
the users simultaneously take part in creating chaos which thus becomes common for 
all the network and therefore recognizable by any participating chaotic generator. At 
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present form, however, this scheme needs in temporal separation of users and other 

hmitations [40,41]. 
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Conclusion 

We hope that at least some parts of the aforesaid material can be useful supplement 
to existing literature on magnetic waves. Anyway, the presented approach, - based 
on first principles only, - well helps to understand and interpret results of numerical 
simulation of linear and non-linear magnetostatic waves and magnetic chaos. This will 
be subject of continuation of this manuscript. 
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